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Abstract 

For thin films of superfluid adsorbed on a disordered substrate, we derive 
within a mean field (Hartree) description of the condensate the equation of 
motion for a vortex in the presence of a random potential. The compressible 
nature of the condensate leads to an effective pinning potential experienced 
by the vortex which is nonlocal, with a long range tail that smoothes out the 
random potential coupling the condensate to the substrate. We interpret this 
nonlocality in terms of images, and relate the effective potential governing 
the dynamics to the pinning energy arising from the expectation value of the 
Hamiltonian with respect to the vortex wavefunction. 
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I. INTRODUCTION 



Pinning of vortices in bulk helium has a long history (see for example Ref. [l[]) but the 
study of vortex pinning in thin films is at a much earlier stage of development: to our knowl- 
edge only one experimental study having been performed — Ref. |2|], although there are some 
aspects addressed in Ref. 0. In this paper we will concentrate on two dimensional pinning 
of idealised (in a sense defined presently) superfluid films which provides two simplifications 
compared to the situation in three dimensions: the ability to change the healing length as 
well as the obvious contrast between point vortices and extended vortex lines. The healing 
length is altered by changing the coverage of the Helium film, although quantifying this 
relationship is not easy. 

Ellis and Li showed that by 'swirling' a gold-plated mylar substrate a remanent 
vorticity could be created, where the density of pinned vortices was 50,000cm~ 2 so that the 
separation (45/i) was of the order of 10 4 larger than the film thickness (3.2nm). Putnam 
et al || have studied the topography of gold films deposited under similar conditions to 
those used by Ellis and Li; they find that the films have a surface of 'rolling hills' with a 
characteristic length scale of 300Ato 3000A. Thus the separation was considerably larger 
than the topographic features which are presumably the pinning agents. In this paper we 
will concentrate on the pinning of single vortices (and their dynamics) which should be 
appropriate under these experimental conditions. However we will confine ourselves to the 
simplest case of monolayer films. 

We wish to derive from first principles the form of the pinning potential which a vortex 
experiences, given the potential in which the helium atoms move due to the substrate. We 
will see that these two quantities are not the same. We will describe the condensate at 
a mean field (Hartree) level using the hydrodynamic representation |J (in terms of the 
two-dimensional density of the film and velocity potential of the flow) and of the resulting 
nonlinear Schrodinger equation (NLSE) which determines the condensate wave function ||. 
We will allow for some effects of compressibility — for instance that the (two-dimensional) 
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density is modified to 'screen' the random potential |7j — and will find an equation of motion 
that bears some relation to the Magnus effect. 

In this very thin film limit, there are two regimes in considering the behaviour of a vortex 
in a random potential. Firstly there is the region, near the centre of the vortex, where the 
largest contributions to the energy density are the kinetic energy of the fluid and the change 
in the films density that is caused by the flow. However at large distances, the dominant 
contribution to the energy density is the response of the density to the random (pinning) 
potential. Thus the problem falls into two parts and the manner in which to match the 
approximate solutions in those regions. The matching is performed by following the analysis 
of Neu || which was constructed to determine the motion of a many vortex system in a 
compressible ideal fluid from first principles, reproducing the Kirchhoff result for widely 
separated vortices. There the two regions are: where single vortex effects are predominant 
(near the centre of each vortex) and where many vortex contributions are important (on 
longer length scales). Other expansions similar in spirit have been for compressible vortex 
rings || and a number of works following Neu, which may be traced from Ref. [[IT]]. 

In the monolayer regime, there is an additional complication: the superfluid-insulator 
transition (see, for example, the recent publications Ref. |TTJ and Ref. pJJ] and references 



therein). However sufficiently far away (a small fraction of a monolayer) from the onset 
of superfluidity, one may describe the helium film as being an 'inert', nonsuperfluid, initial 
layer with a mobile, superfluid, film adsorbed on top The treatment in this paper 

will be of the latter part of the film, where the random potential is the residual one which 
includes any interaction with the inert layer. 

The plan of the paper is as follows: In section 2 we introduce the framework used to 
describe the vortex motion. Using a perturbation calculation, we derive the velocity of 
a vortex under the influence of an external potential, by matching the 'inner' and 'outer' 
solutions in section 3. In the next section, we establish that the dynamics are nonlocal in the 
random substrate potential and explore some of the consequences. In the following section 
we relate the motion to the gradient of the expectation value of the energy. In the final 



section we discuss several issues which emerge from the rest of the paper and conclude. 

II. NOTATION AND FLUID REPRESENTATION OF THE NLSE 

We wish to study the influence of a disordered substrate qualitatively on a superfluid film. 
We therefore pick a very simple interaction between the bosons, namely a point interaction. 
Hence our starting point is the Hamiltonian 

N t2 i N N 

H = T,-^ + ^T, % - + A £ V(T t ) (1) 

i=l i^j i=l 

Here the coordinates and Laplacian are two-dimensional, as we assume motion quantised 
normal to the substrate surface, with all atoms in the lowest state of that motion. A is the 
variance of the random two-dimensional potential V: 

A 2 = ([^(r)] 2 ) (2) 

where the average is over the ensemble of potentials. Note that V itself is dimensionless. 

We will treat the above Hamiltonian within the Hartree approximation, thus we will not 
include effects at the Bogoliubov approximation or beyond: for instance any analogue of 
the roton minimum. The Hartree approximation leads to the time-dependent NLSE which 
governs the condensate motion f)|||: 

-\v 2 <P+\<P\ 2 <P + oV<P = i^ (3) 

Where we are using appropriately scaled units: lengths are measured in units of the healing 
length, £h = (h 2 /Xnm) 1 / 2 ; a = A/n\ is the dimensionless measure of the strength of the 
external potential, and energy is measured in units of nX, which is the Hartree energy or 
chemical potential, n is the average density of particles (equivalent to the condensate density 
in this Hartree case). <fi is normalised to the size of the system, Q: 

I \<f)\ 2 dV = Q (4) 
Jn 
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so that in the absence of a potential, 0=1 everywhere within f2. Note that the speed of 
sound in the condensate, (nA/m) 1 / 2 , is equal to unity with the above choice of units. 

To describe the dynamics of the condensate possessing a vortex structure, it is is natural 
to use the (Hartree) fluid, or Madelung, representation ||; by setting = ^[pe %s , Eqn. ([|) 
is equivalent to the pair of equations: 

|| + V-(pVS) = (5) 

and 

Eqns. (I]) and (§) are respectively, the continuity, and Bernoulli equations describing the 
condensate flow. Density, p, is measured in units of n. 

In the main body of this paper, we will be interested in the behaviour of the flow field 
well beyond the core region, so that the rapid variations in the condensate density and 
velocity field near the core itself are not considered. Within this approximation we can 
neglect the 'quantum pressure' term V 2 v /p/ \fp in the Bernoulli equation ([]), (which we 
justify a posteriori presently) and use the set of equations: 

^ + Vp • VS + P V 2 S = (7) 

and 

^ + -(VSf + p + o-V = (8) 
For this work, an important solution to these equations is that of a vortex in two-dimensions: 



S = 9-t p=l-~ (9) 



where 6 is the polar angle. Note that the term —t in the expression for S is the chemical 
potential, which cancels the term in the asymptotic density in the Bernoulli equation. 
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III. DYNAMICS DUE TO THE SUBSTRATE POTENTIAL 



In the Magnus equation for incompressible point vortices, the presence of a potential 
coupling to the vortex will change its motion from being determined completely by the value 
of the velocity field at the vortex centre. (The latter result, due to Kelvin and Helmholtz, 
has been derived in the context of Landau- Ginzburg theory by Neu ||.) In this section we 
will derive an approximate solution for a vortex in the presence of a potential which couples 
to the helium atoms in the film. This will allow us to derive in the next section the form of 
the potential which applies to the vortex, as against the helium atoms, and hence appears 
in the equation of motion of the vortex. We apply Neu's [§] method of matched asymptotic 
expansions in this section. 

We will perform a linear response calculation using the fluid (Madelung) representation 
of the condensate wave function. This is preferable to using linear response theory directly 
- in terms of the condensate wave function and the Bogoliubov excitations in the presence 
of the vortex. In the latter approach we would need to calculate the matrix elements of, 
for instance, the random potential between the ground state and the excitations. Since the 
wave functions of the excitations around a vortex are only known asymptotically this would 
not be straightforward. 

Before embarking on any calculation, we wish to divide the problem into one of 'inner' 
and 'outer' solutions. In the inner region, the vortex motion is dominant, in terms of both 
density and velocity field, and in the outer region, the random potential is dominant. We 
can estimate the boundary between these two regions to be where the perturbation in the 
density due to flow around the vortex is equal to the perturbation to the density due to the 
random potential; namely at to be a distance r c where: 

— ~ a => r c ~ a- 1 ' 2 (10) 

' c 

and r c is measured in units of the healing length. 

The inner solution can be determined by perturbation theory on Eqns. (]?]) and (^) in 
powers of a, since the random potential is assumed to be weak. 
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For the outer solution, we want to rescale the position and time variables r and t to make 
the dominance of a at large distances manifest. We make the substitutions: 

r^r' = ffr, t ->• t' = a 2 t (11) 

In terms of these new variables, the continuity and Bernoulli equations of (^) and @ become 
(dropping primes): 

^ + Vp-VS + pV 2 S = (12) 

and 

+ ° 2 \(VS) 2 + P + aV = (13) 

The continuity equation remains unchanged but the Bernoulli equation is modified — we may 
neglect the time-dependence of S (apart from the term — t/cr 2 which provides the chemical 
potential, as mentioned below Eqn. (||). In the outer region we expand the density and 
velocity potential as: 

P =1 + <t Pi (t), S = S (r-R) + aS 1 (r) (14) 

Note that the density does not depend on the position of the vortex, as the terms in Eqn. 
(|j"3|) which depend on R are second order in a, similarly we may neglect the time-dependence 
of p as time does not enter Eqn. (|i"3f) to first order in a. In Eqn. (|i"4|) the zeroth order flow 
field retains a vortex flow, so that § V5*o ■ dr = 2n, despite a distortion of the flow as a whole 
due to the linear addition of VSi. 

Substituting the expansion ([14]) into Eqn. (|12"D and (|i"3f) we have to 0(a): 



Pi = -V 

V 2 S l + VS ■ Vpi = (15) 
Solving for Si and denoting the outer solution by S±, we find: 

S>(r) = J G(t - r) 7 ■ W(r') dr' (16) 



where 

G(r-r') = ^log|r-r'| (17) 
Zix 

and we have introduced the more compact notation 7 = VSo(r — R) for the vortex field. It 
is important to note that the variables here are the scaled ones defined in Eqn. ( |TT|) . 

Turning to the inner solution, we revert to the unsealed variables, and solve the original 
equations (0) and (||). We require the solution to be matched to the outer one at distances 
much greater than the core radius, i.e. at dimensionless distances r ^> 1. In the inner region 
we can write to 0(a): 

p = p (r — R) + api , S = S (r - R) + aVSi 

with p given by 

P ° = 1 ~ 2|r-RP 
= l-5po(r) 

(18) 

So Spo is the same as the change in density due to the vortex in the absence of the random 
potential. The fluid equations taken to 0(a) are then: 

Pi = -V + ~/-(R-VS 1 ) 
p V 2 St + 7 • Vpi + Vp • (VSi - R) = (19) 

The eventual aim is to find R in terms of Si; that is, we look for the change in flow caused 
by the potential, which then determines the vortex velocity through advection. Examining 
Eqn. (fLlp , we can inspect the individual terms to see which ones will dominate in the r — > 00 
limit for the matching. Eliminating p x , and omitting terms of order 1/r 3 and higher, we are 
left with 

V 2 Si = 7 • w 
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which can be solved Sf(r) for to obtain 

S<(r) = Jg(t- r') 7 • Vy(r') dr' + 0(1 /r) 

where S x = Sf denotes the inner solution. 

We will now match the inner and outer solutions at a distance r ~ o~~ a , with a > 
chosen to ensure corrections are small. Firstly, we must consider the outer solution back 
in terms of the unsealed variables. But from Eqn. (0), it can be seen that there is in fact 
no change in its form, as lengths only enter S± and the scaling of 7 • V and the element 
of integration dr cancel. Hence we can conclude that matching is perfect up to 0(a). At 
r ~ cr~ a , neglected terms in the equation determining the inner Si are 1/r or a factor of o~ a 
smaller; thus any inconsistency is negligible. We will make a more refined consistency check 
presently. 

Note that (neglecting corrections) Si can be integrated by parts to give: 

Si(r) = - J 7 (r') • VG(r - r') V(r') dr' (20) 

where we assume that the 'surface' integral vanishes (and V ■ 7 = 0). This alternate form 
for Si will be useful in the next section. 

So far the velocity of the vortex, R, has not been determined. To do so, we Taylor 
expand the outer solution fll6|) about r' = R' and check that corrections to this expansion 
are negligible in the matching region. This ensures that the Taylor expansion is consistent 
with the matching ||. Thus we have 

S>(r') = S>(R) + (r' - R) • VS>(R) + 0(r' 2 ) 

In the second term on the right, VSf appears as VSjf is perpendicular to (r' — R'). 

We now appeal to the Helmholtz theorem which states that the vorticity moves with 
the local velocity field (the 'dynamical boundary condition'). Using Eqn. ([TEj) this gives the 
result for R as: 

R = J V R G(R - r) 7 • W(r) dr (21) 
9 



where Vr = d/dH. 

To check consistency, we examine the size of (V — R') • R = o~(r — R) ■ R compared with 
the corrections 1/r and a 2 r 2 (i.e. the 0(r' 2 ) in the Taylor expansion). We require or ^> 1/r 
and <7 2 r 2 for r ~ cr~°. The first condition gives a^ 1 "^ ^> o~ a =>• a > 1/2, while the second 
one gives a^ l ~ a " > ^> cr( 2 ~ 2a ) =>- a < 1. Hence for both sets of corrections to be negligible, we 
need: 

l/a 1/2 < r < 1/(7 
Since this region does exist, the expansion is consistent. 

IV. INTERPRETATION OF VORTEX MOTION IN TERMS OF IMAGES 

We have now deduced the form for R which governs the vortex motion in an external 
potential, and in this section consider some of the consequences and the interpretation of 
Eqn. ([H]). Firstly, note that R is clearly related to the potential in a nonlocal manner. 
Using Eqn. fl2"0|), the expression for R in Eqn. (^TJ) can be rewritten as: 

R = -V-r J ~f(r) ■ VG(R-r)V(r)dr (22) 

We can further manipulate this by using the identity 6 ■ W = (z x W) ■ r where z is 
perpendicular to the plane of motion, so that (|22|) becomes: 

R = -z x J(/j, ■ V R ) V r G(R - r)V(r) dr 

with fi = (R- r)/|R- r| 2 . 
This implies: 

z x R = V R y cfr 

where 

This reproduces the form of the Magnus equation but with an effective potential V e g which 
is nonlocal in the substrate potential with a 1/|R — r| 2 kernel. 
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The most natural interpretation of this result is in terms of images. To indicate the 
relevance of images, note that if the potential were sufficiently large then the density of the 
film would become zero. For simplicity assume that the boundary of the film is a straight 
line. In that case in a standard manner we may take into account the boundary condition 
that the velocity field perpendicular to the film boundary is zero on the boundary by adding 
the velocity field of an image vortex. It is plausible that when the potential causes changes 
in the density, but is not sufficient to drive the density to zero, there will still be features in 
the velocity field which may be attributed to images. Their presence ensures that the flow 
field obeys the continuity equation with terms involving Vp. We will now indicate, in terms 
of a simple example, how Eqn. (|23f) is interpretable in such a manner. 

We now give an example of vortex motion for the case of a repulsive delta function 
potential centred at r a : 

V(r) = 5(r - r a ) (24) 
The Greens function ( |T7D satisfies the identity: 

VVG(r) = {cos29(e x e x - e y e y ) + sm26(e x e y + e y e x )} (25) 



Using Eqns. ( f24|) and (|25|) , the vortex velocity (from (|22|) ) is given by: 

R-r a 



* = T^TTs (26) 



where 6 a is perpendicular to r a — R. 

To interpret this result, let us compare an image approach to the motion of a vortex in 
the presence of an impenetrable circular region centred at r a . Then the magnitude of the 
velocity experienced by the vortex due to the image is of 0(1/ d 3 ) where d is the distance 
between the vortex and its image (this is because the magnitude of the velocity field due 
to the vortex itself at the circle is o(l/d) and hence the dipole moment required to cancel 
out the normal component is of this magnitude, hence the velocity field back at the vortex 
will be 0(1/ d x 1/d 2 ). This has the same functional form as Eqn. (|26[) , leading to the 
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interpretation in terms of a dipole image at the delta function potential, but with a strength 
which is not unity but proportional to the strength of the potential. This has some analogy 
with the difference between images in electrostatics in the case of metals and dielectrics. 
The direction of the motion of the vortex is also consistent: in the image picture the vortex 
would be advected around the circle, in agreement with the angular unit vector in Eqn. 
(P6f). It is easy to check that a line of delta functions produces motion parallel to the line 
(of magnitude 1/d), in analogy with the motion of a vortex in the presence of a wall. 

We note parenthetically that the response of the condensate to a delta function perturba- 
tion without a vortex is of a much simpler form: a delta function change of the condensate. 
This is because the repsonse of the condensate to a perturbation relaxes on a scale of the 
healing length; the fluid description only represents behaviour on scales larger than the 
healing length, hence the perturbation is of a delta function form within the fluid approxim- 
tation. However as we have seen, there is an interesting and nontrivial response on length 
scales in excess of the healing length in the presence of a vortx. 

Having established that the vortex obeys a Magnus equation of motion and hence moves 
parallel to the equipotentials of the effective potential, V e s, we may immediately draw some 
additional conclusions using the work of Trugman and Doniach A consequence of Eqn. 



( P3|) is that the vortex trajectories must be closed, with V e s instead of V as the underlying 
(local) potential. In other words, the vortex travels along the equipotential lines of the 
former. The exception to this is at the saddle points of the effective potential where the 
vortex may 'percolate', allowing in principle an extended orbit across the system. (One 
can think of of 'hills' and 'valleys' representing the potential, and the percolation threshold 
corresponding to 'lakes' which fill up the valleys, to connect). However, the speed of the 
vortex goes to zero (due to the vanishing of the gradient of the potential) at the saddle 
points, so that the vortex transport is in fact pathological. These points have immediate 
parallels with the guiding centre motion of charged particles in the quantum Hall effect 
]TBHi7[, where the existence of percolating paths is related to the mobility edge believed 
to occur in the middle of a Landau sub-band. Obviously the application of these ideas to 
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vortices would only be appropriate when the vortex density was extremely low. 



V. PINNING ENERGY 



So far the discussion has been in terms of the equations of motion. In this section we 
show that we may calculate the pinning energy and moreover the force which occurs in the 
equation of motion of the vortex is the gradient of that potential. 

The expectation value of the energy is: 



Now, substituting in the expansions for the density and velocity potential to 0(a), there 
are four terms which depend on the vortex position contributing to this energy (here Spo(r) 
is defined in Eqn. ([18])): 



where we have taken the vortex core to be situated at the origin. 

We can interpret the separate terms on the right hand side of this equation as follows. 
The first term is just the pinning potential energy associated with the decreased density as 
the core is approached. The second term is the interaction of this decrease in the density 
with the distortion of the condensate due to the substrate potential. The third term comes 
from the change in the kinetic energy due to the substrates effect on the density. The final 
term comes from the change in the kinetic energy due to the distortion of the flow field 
due to the substrate (i.e. the consequence of the vortex flowing around 'obstacles' in the 
substrate). Now the first and second terms cancel as p\ = —V. Thus we are left with the 
final two terms: 



The (inverse) power law nature of 5po means that the pinning is nonlocal as to be expected 
from the discussion in section 4. We may now rewrite the second term using the two- 
dimensional divergence theorem: 



E = \ P (VS) 2 + l -p 2 + aVp 




(27) 




(28) 



13 



Jdrj-VS 1 = J V • (7S1) - J dr V ■ 7S1 
= JdS-jSr- J drV-75i 

= JdS--yS 1 -0 (29) 

as V ■ 7 = 0. 

The surface integral in Eqn. (p9|) will be zero as the velocity field parallel to the boundary 
will vanish due to additional images (which we do not discuss explicitly). The final result is 
thus: 

E = -J -^ 2 Vdr 

i.e. the interaction of the depletion due to the centrifugal forces interacting with the substrate 
potential. 

It should be stressed that the core region, r < 1, will contribute to the pinning energy 
a term of comparable magnitude to the long range part which we are explicitly discussing; 
however the latter yields only a short range correlated potential as against the long range 
correlations which the above has. It is readily checked that minus the gradient of this 
potential is indeed the term on the r.h.s of the equation of motion of the vortex and hence 
is consistent with our considerations of the equation of motion. 

VI. DISCUSSION AND CONCLUSIONS 

In this paper we have only discussed the motion of a single vortex in the presence of a 
random potential. Let us now briefly discuss the relevance of this to the more realistic case 
where there are many vortices present in the system. The vortices may be the remanent 
vorticity which is pinned by the random potential, or it may be that the substrate is rotating 
and hence there is some vorticity with a net sign. 

We concentrate on the case where the disorder is strong compared to the interaction of 
the vortices at the average separation. Of course in the case of a rotating substrate the 
logarithmic nature of the interaction between the vortices will ensure that the density does 
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not fluctuate too much, however the local order in the vortices is determined by the higher 
Fourier components of the random potential, and the shear modulus of the vortex lattice 
which will decrease with density. 

Because of the long-range interactions we must consider the possibility of an analogue of 
a Coulomb gap which occurs for electronic excitations in a highly disordered system JTE|-pO 



The 'gap' refers to the density of states for adding an extra vortex to the system. In the 
case of vortices with a logarithmic interaction, the density of states, n(e), as the energy, e, 
tends to zero, behaves as n(e) ~ e 3 . However the situation with excitations where a vortex 
is moved within the system is more complicated as one is creating a dipole and the density 
of these excitations does not tend to zero at zero energy. The above results all assume that 
the particles (vortices in our case) reside on a lattice; however if the potential is smooth 
then as well as excitations where the vortex is moved between local minima in the overall 
potential (due to the combination of the underlying potential and the effect of the other 
vortices) there are excitations where it remains in the minimum which it started in. It is 
those excitations with which we have been concerned in this paper. 

These excitations would (at least in principle) contribute to the damping of a torsion 
balance. This may be seen readily by considering a vortex in an approximately harmonic 
potential well (V e s(r) — (l/2)Kr 2 ) experiencing a superflow with amplitude A and frequency 
Q in the direction x. The equation of motion for the vortex with position r and unit 
circulation is then: 

(r — A cos QtH) x z = —Kr 

The effect of the superflow is to give a response in terms of the motion of the vortex which is 
resonant if Q = ujq = K. (That is if the period of the superflow is the same as the frequency 
of the vortex moving around the well.) For instance the x-component of the vortex position 
is: 

AM . „ 

x = a cos Uot + a sin u t + — „ sin i It 

\V — Uq 

So those vortices which have associated natural frequencies which are the same as that of 
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the torsion balance will be excited significantly. Of course this frequency is rather low on a 
microscopic scale and so it may be insignificant in practice. 

The extent to which the excitations are single-vortex in their nature is hard to estimate 
in general. However one may start from the assumption of only a single vortex is excited and 
then examine whether the amplitude spreads to other vortices. The single vortex excitation 
has a characteristic frequency involved in precessing around an equipotential of Vgg. Whether 
this may then excite neighbouring vortices depends on whether the natural frequencies of 
the neighbouring vortices are sufficiently similar for a 'resonant' process. If the random 
potential has a sufficiently large variance, as is assumed here, this is very unlikely so most 
excitations will only involve one vortex. Hence the results of this paper will be relevant. 

We have not considered here the effect of the vortices having a nonzero inertial mass. 
The effect of this may be deduced by analogy with the case of a charged particle in a 
magnetic field — which is very similar upon the addition of an inertial mass term. As well as 
the 'guiding centre' motion of the vortex along the equipotentials, there is now some (fast) 
'cyclotron motion'. (A discussion of this effect in classical hydrodynamics is given by Lamb 
). However the value of the inertial mass (and hence the magnitude of the cyclotron 



motion) has been a controversial issue. The suggestion of Baym and Chandler |22[ was that 
it was associated with the virtual mass due to the backflow around a cylinder of radius 
approximately the coherence length. However more recently Duan |23| has argued that the 
mass is larger by a factor of 20-30 due to the finite compressibility of helium. In terms of 
the analogue of the 'cyclotron radius', £ m = (p//i) _1//2 , where fx is the mass of the vortex, 
will only change from being roughly an interatomic distance to being 4-5 times that size. 
In this paper we have nothing to say about corrections to the mean field picture. This 



has been addressed in part by Niu et al. [p^j , where the magnitude of the vortex effective 
mass has been considered from a quantum point of view. The relationship between these 
calculations and the effects due to backflow is not very clear and we will make no more 
comment on these matters. 

In summary we have derived the effective potential that a vortex experiences due to a 
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potential coupled to the underlying bosons. The former is much smoother than the latter, 
with long range tails. We showed that the tails may be interpreted as due to 'images' caused 
by variations in the condensate density. Moreover the effective potential was equal to the 
pinning energy in the expectation value of the Hamiltonian with a Hartree wave-function 
representing the vortex. 
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